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ABSTRACT 

We investigate radial metallicity gradients for a sample of dwarf stars from 
the RAdial Velocity Experiment (RAVE) Data Release 3 (DR3). We select 
a total of approximately 17000 F-type and G-type dwarfs using selection of 
colour, log and uncertainty on derived space motion, calculate for each 
star a probabilistic (kinematic) population assignment to thick or thin disc 
using space motion, and additionally another (dynamical) assignment using 
stellar vertical orbital eccentricity. We additionally subsample by colour, to 
provide samples biased to young thin disc and to older thin disc stars. We 
derive a metallicity gradient as a function of Galactocentric radial distance, 
i.e. d[M / H]/ dR m — —0.051 ± 0.005 dex kpc -1 , for the youngest sample, F- 
type stars with vertical orbital eccentricities ey ^ 0.04. Samples biased to 
older thin disc stars show systematically shallower abundance gradients. 

Key words: Galaxy:abundances - Galaxy: disc - stars: abundances - Galaxy: 
evolution 
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1 INTRODUCTION 



Abundance gradients in galactic discs are an important constraint on the star formation and 
interstellar medium history and the possible secular evolution of stars post-formation in those 
discs. In the Milky Way Galaxy there is exte nsive information establishing a radial gradient 



in young stars and in the inte rstellar medium (IShaver et al 



1983 



2006 



Luck &: Lambert 



Luck. Kovtyukh fc Andrievsky 



201 lh . Values typically derived are d[Fe/H}/dR G = -0.06±0.01 dex 



kpc -1 , within 2-3 kpc of the Sun. Much effort to search for local abundance variations, and 
abundance variat ions with azimuth, show little if any detectable variation in young systems 
( jLuck et al.ll201ll ). limiting the importance of radial gas flows. 

Quantifying the abundance distribution functions and their radial and vertical gradients 
in both the thin and thick discs can be achieved using stellar abund ances, especially those 
from major surveys such as RAdial V elocity Experiment (RAVE, 



Steinmetz et al. 



2006 



Zwitter et al. 



2008 



Siebert et al. 



20111 ) . Several recent analyses of RAVE data have studied 



the metalli city and kinematics of the Galactic discs distribution functions from RAVE data 



( jBurnett et al. 



2011 



Karatas fc Klementl 



2011 



Ruchti et all 



2011 



Wilson et al. 



201lh . 



Even more fundamental is the variation of the abundance gradient with time. Some 
information is available, especially from open star clusters and field stars . The evidence here 



suggests the gradient flattens beyond about 12 kpc Galactocentric radius (ICarnev et al.ll2 005 



Yong. Carney fc Teixera de Almeida 



2005 



2005 



) , as is seen in other spiral galaxies ( iWorthey et al. 



Vlajic. Bland-Hawthorn fc Freeman! 120091 ) . At very outer radii the situation may be 



more complex. 

The observational situation however can be improved. Extant data suggest vertical metal- 
licity gradients in the —0.4 < d[M/H]/dz < —0. 2 dex kpc -1 range for relatively small dis 



tances from the Gal 


actic plane, i.e. 


2003; 


Du et al. 


2004; 


Ak et al. 


2007a) 



1995 



Karaali et al. 



2007af ). For intermediate z distances, where the thick disc is 



dominant, the vertical metallicity gradient is low, d[M/H]/dz = —0. 07 dex kpc 1 and the ra- 
dial g radient is only marginal, —0.02 ^ d[M/H]/dz ^ dex kpc -1 ( IRong. Buser fc Karaali 



20011 ) . There is some evidence that metallicity gradients for relatively short vertical dis- 



tances, z < 2.5 kpc, show systematic fluctuations with Galactic longitude, similar to those 
of the thick-disc scaleheight, which may be interpreted as a common flare effect of the disc 



( 1Ak et al. 



2007b 



Bilir et al. 



2008 



Yaz fc Karaali 



20101 ) . 



These vertical gradients are perhaps a convolution of a time-dependent abundance and 



Radial Metallicity Gradient from RAVE DR3 3 

the well-established stellar age-velocity dispersion relation. Alternatively, they may be more 
consistent with two independent, rather narrow metallicity distributions, one thin disc, one 



thick di sc, with the apparent gradient being 



heights (IGifmore fc Wyse 



1985; 



simp 



Burnett et al. 



201 lh . 



y a reflection of the two different scale 



In this study we analyse stellar abundance gradients from RAVE data, using colour 
bins as a proxy for mean age, to test for time- dependent effects. RAVE is a multi-fibre 
spectroscopic astronomical survey of stars in the Milky Way, which covers just over half of the 
Southern hemisphere, using the 1.2 m UK Schmidt Telescope of the Australian Astronomical 
Observatory (AAO). RAVE's primary aim is to derive the radial velocity of stars from 
the observed spectra for solar neighbourhood stars. Additional information is also derived, 
such as photometric parallax and stellar atmospheric parameters, i.e. effective temperature, 
surface gravity, metallicity and elemental abundance data. This information is important in 
calculating metallicity gradients, which provides data about the formation and evolution of 
the Galaxy. As the data were obtained from solar neighbourhood stars, we have limitations 
to distance and range of metallicity. However, the metallicity measurements are of high 
internal accuracy which is an advantage for our work. A radial metallicity gradient of -0.04 
dex kpc -1 , based on calibrated metallicities from the RAVE DR2 data has already appeared 
in the literature (jKaratas &: Klementl l201ll ). However, this metallicity gradient covers all 
spectral types, and thus a wide age range, which is a disadvantage for a sample of stars with 
distance and metallicity range restrictions. Hence, in this study we will use only the data 
for the F and G spectral type stars. We test for different metallicity gradients for different 
spectral types. 

The structure of this paper is: Data selection is described in Section 2; calculated space 
velocities and orbits of star samples and population analysis are described in Sections 3 and 
4, respectively. Results are given in Section 5 and a summary and conclusion are presented 
in Section 6. 



2 DATA 



The data were selected from the third data release of RAVE (DR3 - ISiebert et al.l 120111 ) . 
RAVE DR3 reports 83072 radial velocity measurements for stars with 9 ^ / ^ 12. This 
release also pro vides stellar atmosp heric parameters for 41672 spectra representing 39833 



individual stars ( ISiebert et al.ll201ll ). The accuracy of the radial velocities is high, marginally 
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improved with DR3: the distribution of internal errors in the radial velocities has a mode of 
0.8 km s _1 and a median of 1.2 km s" 1 , while 95 per cent of the sample has an internal error 
smaller than 5 km s -1 . The uncertainties for the stellar atmospheric parameters are: 250 K 
for effective temperature T e ff, 0.43 dex for surface gravity logg and 0.2 dex for [M/H\. While 
RAVE supports a variety of chemical abundance scales, we use here just the public DR3 
values. Since anticipated gradients are small, this provides a well-defined set of parameters 
for analysis. 

A surface gravity constraint 4 < logg ^ 5 was applied to the 83072 star sample to 



obtain a homogeneous dwarf-star subsample with accurate data. Distances o 



stars were obtained using the absolute magnitude calibration of 



Bilir et al 



the s ubsample 



(120081 ). whereas 



the reddening values were obtained iteratively, f ollowing published met hodology (for more 



detailed information regarding the iterations see ICoskunoglu et al. 
therein) . 



2011 



and the references 



Our aim here is to take advantage of RAVE's focus on high Galactic latitude to isolate, 
statistically, samples dominated by thin disc and by thick dis c stars. Thus, after dereddenin g, 
stars were split into F and G spectral type subsamples using IStraizys &: Lazauskaitd (120091 ) 's 
criteria, i.e. 0.09 < (J-H) < 0.29 and 0.29 < (J - H) ^ 0.39 leaving us with 17768 stars. 
Out of these 17768 stars 10661 are F- and 7107 are G-type dwarfs, while the sample median 
distances are 326 and 272 pc, respectively. 



3 SPACE VELOCITIES AND ORBITS 



We combined the distances estimated in Section 2 with RAVE kinematics and the avail- 
ab le proper motions, app l ying the (standard) algorithms and the transformation matrices 
of I Johnson fc Soderblom ( 1987 ) to obtain their Galactic space velocity components (U, V, 
W). In the calculations epoch J2000 was adopted as described in the Internatio nal Celestial 
Reference System (ICRS) of the Hipparcos and Tycho-2 Catalogues (jESAlll997l ). The trans- 
formation matrices use the notation of a right handed system. Hence, U, V and W are the 
components of a velocity vector of a star with respect to the Sun, where U is positive to- 
wards the Galactic centre (/ = 0°, b = 0°), V is positive in the direction of Galactic rotation 
(/ = 90°, b = 0°) and W is positive towards the North Galactic Pole (b = 90°). 

Correction for differential Galactic rotation is necessary for accurate determination of 
U, V and W velocity components. The effect is proportional to the projection of the dis- 
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tance to the stars onto the Galac tic plane, i.e. the 



by Galactic differential rotation ( IMihalas fc Binney 



W ve. 



ocity component is not affected 



Mihalas &: Binney 



19811) . We applied the procedure of 



( 1l98ll ) to the distribution of the sample stars and estimated the first or- 
der Galactic differential rotation corrections for U and V velocity components of the sample 
stars. The range of these corrections is —25.71 < dU < 16.92 and —1.57 < dV < 2.29 km 
s _1 for U and V, respectively. As expected, U is affected more than the V component. Also, 
the high values for the U component show that corrections for differential Galactic rotation 
can not be ignored. 

The uncertainty of the space velocity components U err , V err and W err were computed by 
propagating the uncertainties o f the proper motions, d i stanc es and radial velocities, again 
using a (standard) algorithm by I Johnson & Soderbloml (119871 ). Then, the error for the total 



space motion of a star follows from the equation: 
s 2 =U 2 +V 2 + W' 



r2 
err' 



(1) 



The median and standard deviation for space velocity errors are S err = 8.01 km s" 1 and 
s = 16.99 km s _1 , respectively. We now remove the most discrepant data from the analysis, 
knowing that outliers in a survey such as this will preferentially include stars which are 
systematically mis-analysed binaries, etc. Astrophysical parameters for such stars are also 
likely to be relatively unreliable. Thus, we omit stars with errors that deviate by more than 
the sum of the standard error and the standard deviation, i.e. S err >25 km s _1 . This removes 
760 stars, ~4.3 per cent of the sample. Thus, our sample was reduced to 17008 stars, those 
with more robust space velocity components. After applying this constraint, the median 
values and the standard deviations for the velocity components were reduced to (U erT , V err , 
W err ) = (3.93 ± 3.10, 3.59 ± 2.79, 3.16 ± 2.75) km s -1 . 

To complement the chemical abundance data, accurate kinematic data have been ob- 
tained and used to calculate individual Galactic orbital parameters for all stars. In order to 



calculate t 
literature ( 


lose parameters we used standard g 


ravitational potentials well-described 


m the 


Miyamoto & Naeai 


1975; 


Hernauist 


1990; 


Johnston, Spereel & Hernauist 


1995; 


Dinescu. Girard & van Altena 


19991) 


to estimate orbital elements of each of the sample stars. 



The orbital elements for a star used in our work are the mean of the corresponding orbital 
elements calculated over 15 orbital periods of that specific star. The orbital integration 
typically corresponds to 3 Gyr, and is sufficient to evaluate the orbital elements of solar 
neighbourhood stars, most of which have orbital periods below 250 Myr. 
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Solar neighbourhood velocity space inclu des well-estab 
ble c l assic moving groups or stellar streams 



199 



0: 



Nordstrom et al. 



Dehnen 



1998 



ished substructures that resem- 



Sku lian. Hearnshaw fc Cottrell 



200J). iFamaev et all ( 120051 ): iFamaev. Siebert fc Jorissenl (|2008j) and 



Pompeia et al.l (1201 If ) show that, although these streams include clusters, after which they 
are named, and evaporated remnants from these clusters, the majority of stars in these 
streams are not coeval but include stars of different ages, not necessarily born in the same 
place nor at the same time. They argue these streams are dynamic al (resonant) in origin, 



proba bly related to dynamical perturbations by transient spiral waves (iDe Simone. Wu fc Tremaine 



20041 ). which radially migrate stars to specific regions of the UV-pl&ne. Stars in a dynamical 



stream just share a common velocity vector at this particular epoch. These authors further 
point out the obvious and important point that dynamical streams are kinematically young 
and so integrating backwards in a smooth stationary axisymmetric potential the orbits of 
the stars belonging to these streams is non-physical. Does this fundamentally invalidate our 
calculations? 



Famaey et al.l (120051 ) assigned probabilities for each star in their sample to belong to 



different kinematic groups. They found that stars not belonging to dynamical streams (young 



giants and the smoot h background) make up the majority (7 per cent) o 



Seabroke et al 



(120081 ) 's fig. 10 illustrates that the majority of IFamaev et al. 



;heir sample . 



([gOOS^ s stars 



are within ±250 pc of the Galactic plane, with the distribution centred on the Galactic 
plane. Our RAVE sample o f F-G dwarfs are at similar line-of-sig ht distances f r om th e Sun 

(120 Ill y s fig. 



as the 



Famaey et al. 



Siebert et al. 



(120051 ) 's K-M giants (200-400 pc). However 
15 shows that RAVE stars are selected to avoid the Galactic plane (|6| > 10° ) . This means 
our R AVE sample will include stars further from the Galactic plane than the IFamaev et al. 



(120051 ) sample. Dynamical perturbations by transient spiral waves are strongest closest to the 
Galactic plane so there will be fewer dynamical stream stars in our RAVE sa mple. While we 



could in principle assign our stars to different kinematic groups, compared to IFamaev et al. 



(120051 ) 's. it is probable that our sample has more than the local 70 per cent of stars whose 
orbital parameters can adequately be deter mined from the static Milky Way disc potential 



and not influenced by transient spiral waves. 



Famaey et al 



(120051 ) 's space velocities are also 



much more accurate than our RAVE values, so considering all factors, we have not attempted 
to remove dynamical stream stars from our sample. Contamination is unlikely to affect our 
statistical results. 

To determine a possible orbit, we first perform test-particle integration in a Milky Way 
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potential which consists of a logarithmic halo of the form 



$haio(r) =ugln 1 + 



d? 



(2) 



with vq = 186 km s 1 and d — 12 kpc. The disc is represented by a Miyamoto-Nagai 
potential: 

GM d 



$disc(-R,2) 



(3) 



with M d = 10 11 M 0) = 6.5 kpc and fe^ = 0.26 kpc. Finally, the bulge is modelled as a 
Hernquist potential 



bulge 



(r) 



GM h 



(4) 



r + c 

using Mb = 3.4 x 10 10 M & and c = 0.7 kpc. The superposition of these components gives 
quite a good representation of the Milky Way. The circular speed at the solar radius is ~ 220 
km s _1 . Plsr — 2.18 x 10 8 years is the orbital period of the LSR and V c = 222.5 km s _1 
denotes the circular rotational velocity at the solar Galactocentric distance, R® = 8 kpc. 

For our analysis of gradients, we are interested in the mean radia l Galactocentr i c dist ance 
(R m ) as a function of the stellar population and the orbital shape. I Wilson et al.l (120111 ) has 
analysed the radial orbital eccentricities of a RAVE sample of thick disc stars, to test thick 
disc formation models. Here we focus on possible local gradients, so instead consider the 
vertical orbital eccentricity, ey. R m is defined as the arithmetic mean of the final perigalactic 
(Rp) and apogalactic (R a ) distances, whereas ey is defined as follows: 



( | Z max | -\- | Z m i n | ) 

Rm, 



(5) 



where R m = (R a + R p )/2 (jPaulil 120051 ) . Due to z-excursions R p and R a can vary, however 



this variation is not more than 5 per cent. 



4 POPULATION ANALYSIS 



4.1 Classification using space motions 



Bensby et al. 



(120051 ) was used to 



The procedure of iBensby. Feltzing fc Lundstroml (120031 ): 
separate sample stars into different populations. This kinematic methodology assumes that 
Galactic space velocities for the thin disc (D), thick disc (TD), and stellar halo (H) with 
respect to the LSR have Gaussian distributions as follows: 
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f(U, V, W) = k x exp^ 

where 

fc - 1 

lsr^Vlsr^W 



TI 2 

U LSR 



LSR 



V ) 2 

* asym ) 



2a 



Ur 



2a, 



WISH ' 

2CT W L SR, 



(6) 



(7) 



normalizes the expression. For consistency with other analyses au LSR , ay LSR and aw LSR were 
adopted as the characteristic velocity dispersions: 35, 20 and 16 km s _1 for thin disc (D); 



67, 38 and 35 


cm s 


(Bensbv et al. 


2003) 



' asym 



is the asymmetric drift: -15, -46 and - 220 km s 1 for thin disc , 



thick disc and halo, respectively. LSR velocities were taken from ICoskunoglu et al.l (120111 ) 
and these values are (U, V, W) LSR = (8.83 ± 0.24, 14.19 ± 0.34, 6.57 ± 0.21) km s _1 . 

The probability of a star of being "a member" of a given population is defined as the 
ratio of the f(U, V, W) distribution functions times the ratio of the local space densities for 
two populations. Thus, 



TD/D 



X. 



X 



I'D fm 
X 



D 



f. 



TD/H 



X. 



D 



td frp 

X , 



(8) 



are the probabilities for a star being classified as a thick disc star relative to it being a 
thin disc star, and relative to it being a halo star, respectively. Xp, Xtd and X^ are the 
l ocal space densities for thin disc, thick disc and halo, i.e. 0.94, 0.06, and 0.001 5, respectively 



Robin et al. 



1996 



Buser. Rong fc Karaalilll999l ) . We followed the argument of 



Bensbv et al. 



(120051 ) and separated the sample stars into four categories: TD/D ^0.1 (high probability 
thin disc stars), 0.1 < TD/D ^ 1 (low probability thin disc stars), 1 < TD/D 10 (low 
probability thick disc stars) and TD/D > 10 (high probability thick disc stars). Fig. 1 
shows the U — V and W — V d i agram s as a function of spectral types and population types 



defined by using 



Bensbv et al. 



(120031 ) 's criteria. It is evident from Fig. 1 that the kinematic 
population assignments are strongly affected by space-motion uncertainties. 

Using Fig. 1, 15270 and 1142 stars of the sample were classified as high and low proba- 
bility thin disc stars, respectively, whereas 287 and 309 stars are low and high probability 
thick disc stars (Table 1). The relative number of high probability thick disc and thin disc 
stars (2 per cent) is evidently very much lower than the number expected, especially for a 
sample biased to intermediate Galactic latitudes. 
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Figure 1. U — V and W — V diagrams of F- and G-type stars applying iBensbv et al,l j 20031) 's population classification criteria. 
It is apparent that space motion uncertainties remain significant, even for this nearby sample. 



Table 1. Population frequency distribution of sample stars according to Bensby et al. (2003) 's kinematical criteria. TD/D 
ranges are explained in the text. 



Number of stars 


Spectral Type 


TD/D < 0.1 


0.1 < TD/D < 1 


1 < TD/D < 10 


TD/D > 10 


Total 


F & G type stars 


15270 


1142 


287 


309 


17008 


F type stars 


9566 


484 


97 


121 


10271 


G type stars 


5704 


658 


190 


185 


6737 



4.2 Classification using stellar vertical orbital shape 

Both radial and vertical orbital eccentricities contain valuable information: here we consider 
the vertical orbit shape. Vertical orbital eccentricities were calculated, as described above, 
from numerically-integrated orbits. We term this the dynamical method. Fig. 2 shows that 
the distribution function of ty for F- and G-type stars, considered separately or in total, 




0.05 



0.1 



0.15 



Figure 2. Three vertical eccentricity distributions: (a) for F- and G- type stars, (b) for F-type stars and (c) for G-type stars. 



are not consistent with being a single Gaussian distribution. A two-Gaussian model however 
does provide an acceptable fit. Hence, we separated our sample into three categories, i.e. 
stars with e v ^ 0.04, 0.04 < e v ^ 0.1 and e v > 0.1, for each spectral type and fitted their 
metallicities to their mean radial distances (R m ) in order to investigate the presence of a 
metallicity gradient for RAVE stars. We provide Table 2 electronically, which includes stellar 
parameters from RAVE DR3, calculated kinematical and dynamical parameters and stellar 
population for the entire sample. The categories and the number of stars corresponding to 
those categories are given in Table 3. 



Table 2. Stellar atmospheric parameters, astrometric, kinematic and dynamic data for the whole sample: (1): ID, (2): Designation, (3-4): Equatorial coordinates in degrees, (5): T e ff in 
K, (6) log <?/ (cms 2 ) , (7): Calibrated metallicity [M/H] (dcx), (8-9): Proper motion components in mas/yr, (10): d in kpc, (11): Heliocentric radial velocity in km s _1 , (12-17): Galactic 
space velocity components, and their respective errors in km s _1 , (18-19): Pcrigalactic and apogalactic distance in kpc, (20-21): Minimum and maximum distance from the Galactic 
plane in kpc, (22): TD/D ratio as mentioned in the text. 
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0) 
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-094053 
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077375 


-9.681528 
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0.10 


-1.60 


-7.50 


0. 
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-5.4 


13.34 


3.40 


3.03 


3.32 


8.27 


1 


.75 


7 


678 


8.475 


-0. 


351 


0. 


351 





.006 


2 


.T000020 


3- 


■ 520757 


0. 


084458 


-52 


132528 


5761 


4. 


06 


0.06 


-12.00 


5.10 


0. 


289 


30.3 


33.37 


6.13 


18.20 


6.11 


-20.21 


2 


.97 


7 


.474 


9.708 


-0. 


386 


0. 


394 





018 
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.7000022 
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-265716 


0. 


094250 


-26 


954694 


6080 


4. 
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0.49 


-6.10 
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0. 
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12.6 


19.15 
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-3.00 


3 
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0. 
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0. 
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.747 
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473 
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474 
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.7000026 


3- 


-394148 
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19 
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0. 
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Table 3. Vertical eccentricity frequency distribution of sample stars. 



Number of stars 


Spectral Type 


e v s: 0.04 


0.04 < e v < 0.1 e v 


> 0.1 


Total 


F type stars 


6923 


3145 


203 


10271 


G type stars 


4223 


2264 


250 


6737 
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Figure 3. Metallicity distributions for: (a) F- and (b) G-type stars. 

5 RESULTS 

5.1 Metallicity Gradients from the Kinematical Population Assignment 
Method 

The metallicity distribution functions of our final sample, divided into two spectral types, 
are shown in Fig. 3. We note these are RAVE DR3 metallicities [M/H] and not standard 
[Fe/H] values. As seen from the figure, the metallicity distribution extends from -0.7 to 
+1 dex. We now consider the metallicities as a function of the mean orbital Galactocentric 
radial distance (R m ) for each different population, i.e. for F and G-type stars and test for the 
presence of a metallicity gradient for each population. We fitted the distributions to linear 
equations, whose gradient is any metallicity gradient, d[M/H]/dR m . The results are shown 
in Fig. 4. The metallicity gradients in all panels of Fig. 4 are either small or consistent with 
zero. The best determined values are for high probability thin disc stars, where the gradient 
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Table 4. Radial metallicity gradients for F- and G-type stars evaluated from kinematical and dynamical data. Ranges of TD/D 
and e v are explained in the text. 



Population type 


F type stars 
d[M/H]/dRm 
(dex kpc - ) 


Sample size 


G type stars 
d[M/^]/dR m 
(dex kpc ) 


Sample 


TD/D ^ 0.1 


-0.043±0.005 


9566 


-0.033±0.007 


5704 


0.1 < TD/D < 1 


-0.024±0.007 


484 


-0.007±0.009 


658 


1 < TD/D < 10 


-0.042±0.017 


97 


O.OOOiO.010 


190 


TD/D > 10 


0.016±0.011 


124 


0.010±0.009 


185 


e v s: 0.04 


-0.051±0.005 


6923 


-0.020±0.006 


4223 


0.04 < e v s: 0.1 


-0.020±0.005 


3145 


-0.004±0.005 


2264 


e v > 0.1 


0.016±0.012 


203 


0.037±0.016 


250 



is d[M/H}/dR m is -0.043 ±0.005 dex kpc" 1 for F-type stars, and -0.033 ±0.007 dex kpc" 1 
for G-type stars (Table 4). Interestingly, the metallicity gradient for high probability thick 
disc stars is consistent with zero. However, we treat this value with caution, for the number 



of stars are only 124 and 185 for F and G-types, resp ectively, for t 
1 reminds us that space motion errors are important, 
thick disc stars in more detail. 



Wilson et al. 



l is pop ulation, and Fig. 
(hoij l discusses RAVE 



5.2 Metallicity Gradients from the Dynamical Population Assignment 
Method 

We divide the data into three bins of vertical orbital eccentricity, and tested for any depen- 
dence of the metallicity versus mean radial Galactocentric distance, R m , again using linear 
fits. The results are presented in Fig. 5 and Table 4. We detect a significant metallicity 
gradient for high-probability F-type thin disc stars of d[M/H]/dR m = —0.051 ± 0.005 dex 
kpc -1 . For G-type stars the metallicity gradients are rather lower being consistent with zero 
(Fig. 5 and Table 4). 

These results are consistent with there being somewhat steeper abundance gradients in 
younger stars than in older. We test that further in the next section. 



5.3 Metallicity Limitation for Metal Rich Stars 

The RAVE pipeline derives metallicity, as any other parameter, using a penalized \ 2 tech- 
nique by finding an optimal match between the observed spectrum and a spectrum con- 
structed from a library of pre-computed synthetic spectra. The metallicity results match if 
a simila r analysis method i s used. The results of the analysis using an independent \ 2 pro- 
cedure ( IMunari et al.l 120051 ) yield metallicities which are entirely consistent with the RAVE 
pipeline results, i.e. mean offset of 0.04 ± 0.02 dex and a standard deviation of 0.17 dex. 
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Figure 4. R m — [M/H] diagrams for sample stars as a function of spectral types and populations. 



RAVE metallicities as derived from the RAVE pipeline are part of a self consistent native 
RAVE system of stellar parameters which are tied to a x 2 analysis using a library of Kurucz 
template spectra. 

However, other spectral methods which derive metallicities from the strengths of indi- 
vidual spectral lines and not from a \ 2 match of synthetic and observed spectra, do not 
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Figure 5. R m — [M/H] diagrams for sample stars as a function of spectral type and vertical orbital eccentricity. 



yield very consistent results with those of the RAVE pipeline. As usual, calibration of an 
internal abundance scale onto standard external systems requires special consideration. One 
can see the following trends in a comparison of metallicities derived by other methods and 
the RAVE pipeline: i) The difference between RAVE and reference metallicity increases with 
an increased a-enhancement, in the sense that RAVE values become too metal poor, ii) The 
difference is larger at lower metallicities. iii) The difference is larger for giants than for main 
sequence stars, though the variation is much weaker than a-enhancement or metallicity. iv) 
The d ifference does not de pend on temperature. 

(j2008l ). the RAVE metallicities ([m/H]) derived by the x 2 method were 



In 



Zwitter et al. 



calibrated to the metallicities ([M/H]) which are in line with the metallicity system of the 
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datase ts in situ, such as those appearing in the Asiago Observatory and 



Soubiran fc Girard 



(120051 ) catalogues. The final calibration is as follows: 



[M/H] = 0.938[m/#] + 0.767[a/Fe] - 0.064 log g + 0.404. (9) 
The classical indicator for the metal abundance is the iron abundance, [Fe/H]. The following 



relatio n between the calibrated metallicity and the iron abundance taken from 



Zwitter et al. 



( 120081 ) produces a limitation for metal rich stars: 



[M/H] = [Fe/H] + 0.11 [1 ± (1 - e -3.6|[Fe/fl]+0.55|^ (10) 

where the plus sign applies for [Fe/H] < —0.55 dex and the minus sign otherwise. As the 
observed upper limit of metal rich stars is about [Fe/H] = +0.5 dex, this is the case for 
calibrated metallicity as well, i.e. approximately [M/H] ^ 0.5. 



5.4 Metallicity Gradients from Blue Stars with [M/H] ^ +0.5 

Our analyses above (Table 4) show that there are differences between the metallicity gra- 
dients estimated from F and G spectral type stars, in the sense that abundance gradients 
derived from RAVE F type-stars are steeper than the metallicity gradients from G type 
stars. We test this result further by separating the data into two spectral-type subsamples, 
F0-F3 and F4-G9, and estimate metallicity gradients for these subsamples. The results are 
shown in Fig. 6. The significant result is that the metallicity gradients for high probability 
thin disc stars (TD/D ^0.1 in Table 4) became steeper, i.e. —0.054 ± 0.008 dex kpc -1 and 
—0.050 ± 0.015 dex kpc -1 for stars of population type F0-F3 and F4-G9 respectively. 

As noted in Section 5.3 the RAVE [M/H] distribution does extend to very high abun- 
dance values, compared to [Fe/H]. In order to test if the stars with most extreme [M/H] 
values have systematically unreliable astrophysical parameters,we excluded stars with RAVE 
[M/H] ^ +0.5 dex from the sample. We then re-estimated metallicity gradients for the two 
subsamples cited above, F0-F3 and F4-G9. The result, given in Fig. 7 confirms our apparent 
trend, in the sense that the metallicity gradient estimated from the bluer stars, —0.065+0.018 
dex kpc -1 , is steeper than that obtained from the later spectral type stars, —0.025 ± 0.008 
dex kpc -1 . 
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Figure 6. R m -[M/H] diagrams for two subsamplcs. (a) for F0-F3 and (b) F4-G9 spectral types. 

6 DISCUSSION AND CONCLUSION 

We have used the RAVE DR3 to identify stars classified as dwarfs, further excluding cool 
stars, and those with the most uncertain space motions. We have then obtained a homoge- 
neous sample of dwarfs defined by: 1) 4 < logg ^ 5, 2) 0.09 < (J — H) ^ 0.39, and 3) total 
space velocity error S err ^ 25 km s -1 . 10271 stars of this sample are F-type stars, whereas 
6737 of them are G-type stars. For each star we calculated distances, total space motions, 
and integrated the stellar orbit in a Galactic potential to derive both mean Galactocentric 
distance and stellar orbital shape. We used the calibrated RAVE DR3 metallicities and these 
mean radial distances to investigate the presence of a radial metallicity gradient, dividing the 
sample into a variety of subsamples. These subsamples were defined using probabilistic pop- 
ulation assignment, considering in turn space motions and stellar Galactic orbital properties 
- what we term "kinematic" and "dynamic" methods. Kinematic properties allowed stars to 
be separated into four populations, i.e. high probability thin disc stars, low probability thin 
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Figure 7. R m -[M/H] diagrams for the same subsamples in Fig. 6, but with mctallicity restriction with [M/H] < +0.5 dex. 
This selection leads to a steeper mctallicity gradient for early spectral type stars (panel a) and a flatter gradient for later 
spectral type stars (panel b). 

disc stars, low probability thick disc stars and high probability thick disc stars, although 
only the high probability thin disc category retained sufficiently large sample size for more 
detailed consideration. In the dynamically-defined sample, stars were separated into three 
subsamples according to their vertical orbital eccentricity, i.e. e v ^ 0.04, 0.04 < e v < 0.10 
and e v > 0.10. 

In all cases, we derive significant metallicity gradients for the samples which are - sta- 
tistically - dominated by the youngest thin disc stars. We derive significant and marginally 
shallower gradients for samples which are - statistically - dominated by somewhat older thin 
disc stars. We do not detect any gradient in the very small samples of stars categorised as 
thick disc. The accuracy and resulting sample size of our thick disc population assignments 
is too poor for any robust statements to be made. 

Our radial metallicity gradients for high-probability thin disc stars, defined using either 
of our two population classification methods, (TD/D ^0.1 and e v ^ 0.04 in Table 4) can 
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be compared to values in the literature. Our F star gradient is —0.043 ± 0.005 dex kpc -1 . 
This is somewhat less than the d[F e/H]/dRc, = —0.06 ± O.Oldex kpc -1 obtained for young 



stars (Cepheids 
gradie nt found by 



Luck et al 



2011 



), but is c onsistent with the thin di sc radial metallicity 



Karatas fc Klementl (120111 ) using earlier RAVE data. IKaratas fc Klement 



(120111 ) considered RAVE DR2 data, combining all spectral types, and using the RAVE 
"calibrated" metallicity scale, to derive d[M/ H]/dR = —0.04 dex kpc -1 . Any difference 
between the RAVE and Cepheid gradients may be due to calibration differences. Our thin- 
disc G star sample, defined using our kinematic approach, TD/ D ^ 0.1 produces a g r adient 
(—0.033 ± 0.005 dex kpc -1 ) which is also consistent with the IKaratas fc Klementl (120111 ) 
value. 

When we define our F and G star samples using our orbital shape (dynamic) definition 



e v ^ 0.04 the resulting abundance gradients 



(-0.0 51 ± 0.005 and -0.020 ± 0.006 dex kpc 



agree less well with the IKaratas fc Klementl (1201 if ) value. Interestingly they are closer to the 
value derived from RAVE DR2 data using RAVE DR2 calibrated [M/H] (-0.04 dex kpc -1 ) 
than the uncalibrated [M/H] (—0.07 dex kpc -1 ). The difference between the values derived 
from RAVE calibrated [M/H] is partly due to improved statistics, RAVE DR3 having wl.5 
times more stars than does DR2, but also warns that RAVE metallicities remain to be 
robustly calibrated onto an [Fe/H] metallicity scale. 



The radial iron gradient for thin disc stars (4 < age < 6 Gyr) found by 



Nordstrom et al 



Schonrich fc Binney 



(2004) is —0.099 ± 0.011 dex kpc -1 . The chemical evolution model of 
( 120091 ) was tuned to provide an excellent fit to this gradient. In the model the gradient 
is caused by radial migration of stars and flow of gas through the disc. This gradient is 
at least a factor of two steeper than our values, presented in Table 3. However, as dis- 
cussed in Section 3, our sample of RAVE stars are further from the Galactic plane than the 
(120041 ) stars, so vertical abundance gradients may - or may not - be an 



Nordstrom et al. 



additional parameter. At face value t he amplitude of time-dependent effects in models such 



as that of 



Schonrich &: Binneyl (120091) may need to be re -tuned to have a smaller effect for 

(120041 ) 's sample because radial migration 



stars further from the plane than 



Nordstrom et al. 



of stars and flow of gas through the disc appears to be less effective further from the plane. 

The steepest metallicity gradient which we measure, i.e. d[M/H]/dR m = —0.051 ±0.005 
dex kpc -1 , corresponds to the subsample of F-type stars defined using orbital shape, with 
eccentricities ey ^ 0.04. The number of stars in this interval is 6923 (67 per cent of the 
whole sample). It is interesting to note that this value is not only consistent with Cepheid 



Rahimi et al. 



( 1201 lh from 



20 Co§kunoglu et al. 

values, but agrees with the value (—0.050 dex kpc" 1 ) predicted by 
their cosmologically simulated galaxy for 'intermediate' disc stars (7 < age < 10 Gyr). In 
their model negative radial metallicity gradients are due to inside-out formation of the disc. 
Our stellar sample is however biased to very much younger ages than those modelled by 
Rahimi et all ( 120111 ). 

Although our samples are small and uncertainties large, we do not detect any significant 
abundance gradient in thick disc stars, using either definition. Our thick disc radial best- 
fit metallicity gradients (TD/D > 10 and e v > 0.1 in Table 4) are however (marginally 



significantly) shallower than our gradients 
thick disc stars (age > 10 Gyr) found by 



or thin disc star s. The radial iron gradient for 



Nordstrom et al. 



feOQJ) is +0.028 ± 0.036 dex 



kpc . Our values are closer to the +0.01 ±0.04 dex kpc found by 



Ruchti et al. 



(1201 lh for 



the metal-poor thick disc. Interes t ingly t hese values agree with chemical ey o 



Chiappini. Matteucci fc Grattonl ( 119971 ); IChiappini. Matteucci fc Romano! (120011 ) (+0.01 — 



ution models of 



0.03 dex kpc 1 ). These models decompose the disc into radial annuli that neither exchange 
gas nor stars. Thick disc stars are, on average, further from the Galactic plane than thin disc 



stars and so should be relatively independent of radial mixing (unless radial mixing actua 



ly 



evolves a thin disc star into a thick disc star as suggested by lSchonrich fc Binneyl (120091 ) ). 
However, within our large errors, thick disc stars show no radial gradient at all. This trend 
is different to what is seen in the thin disc, suggesting these stars do not share the same 
origin as the thin disc. 

Metallicity gradients we have estimated from F-type stars are steeper than the ones 
obtained from G-type stars, for a given population (Table 4). Given our sample size, the only 
robust metallicity gradients we determine are —0.043 + 0.005 dex kpc" 1 and —0.051 + 0.005 
dex kpc" 1 respectively, for stars classified as high probability thin disc, i.e. for population 
types labeled with TD/D ^ 0.1 and ey ^ 0.04. Separation of the sample into two colour 
sub-samples, F0-F3 and F4-G9 here, provides a steeper metallicity gradient, —0.065 + 0.015 
dex kpc" 1 , for the earlier spectral type stars. Other subsamples from the RAVE DR3 dwarf 
sample are absolutely small in number, and have relatively large uncertainties. The RAVE 
DR3 dwarf star sample probes stars with orbits with mean Galactocentric radii within 3 kpc 
of the Solar Galactocentric radius. The stars with spectral types F0-F3 which provide the 
steeper metallicity gradient of —0.065 ± 0.015 dex kpc -1 , make up only 3 per cent of the F 
and G type stars with [M/H] < 0.5 dex. 

From our analysis, we may conclude that the RAVE DR3 data may be described as two 
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different subsamples, i.e. a thin disc sample biased to young ages, with a detected metallicity 
gradient, and a thin disc sample biased to somewhat older ages, in which we do not detect 
any metallicity gradient. 
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